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Abstract 

Electron captures on nuclei play an important role in the dynamics of the collapsing 
core of a massive star that leads to a supernova explosion. Recent calculations of 
these capture rates were based on microscopic models which account for relevant 
degrees of freedom. Due to computational restrictions such calculations were limited 
to a modest number of nuclei, mainly in the mass range A = 45-110. Recent super- 
nova simulations show that this pool of nuclei, however, omits the very neutron-rich 
and heavy nuclei which dominate the nuclear composition during the last phase 
of the collapse before neutrino trapping. Assuming that the composition is given 
by Nuclear Statistical Equilibrium we present here electron capture rates for col- 
lapse conditions derived from individual rates for roughly 2700 individual nuclei. 
For those nuclei which dominate in the early stage of the collapse, the individual 
rates are derived within the framework of microscopic models, while for the nuclei 
which dominate at high densities we have derived the rates based on the Random 
Phase Approximation with a global parametrization of the single particle occupa- 
tion numbers. In addition, we have improved previous rate evaluations by properly 
including screening corrections to the reaction rates into account. 
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1 Introduction 



Near the end of their lives, the core of a massive star (M > IOMq) consists 
predominantly of iron and its nuclear neighbors, the "iron group" elements. 
Since nuclear burning processes continue in the layers above, this iron core 
grows, approaching the Chandrasekhar mass limit, and begins contracting at 
an ever increasing rate. In the inner regions of the core, this collapse is subsonic 
and homologous, while the outer regions collapse supersonically. With densi- 
ties larger than 10^ g/cm^, both transport of electromagnetic radiation and 
heat conduction are very slow compared with the time-scale of the collapse. 
The majority of the energy is transported away from the core by neutrinos, 
originating mostly from electron captures on protons and nuclei pQ. Escaping 
neutrinos reduce the entropy of the core, forcing nucleons to remain as part 
of nuclei, until the neutrinos become trapped at densities 10^^-10^^ g/cm^. 
/3~-decay processes are effectively blocked by the electron degeneracy, causing 
nuclei in the core to become progressively more neutron-rich. 



Once the density in the innermost part of the collapsing star exceeds that of 
nuclear matter, the collapse is halted by the short-range repulsion of the nu- 
clear interaction. The supersonically infalling matter of the outer core bounces 
off this extremely stiff inner part, reversing its velocity and forming a shock 
wave that propagates outwards through layers of lower density. The shock 
weakens as it loses energy to nuclear dissociation of the shocked matter and 
once the shock reaches the neutrino emission surface (or neutrinosphere), the 
remaining energy of the shock is carried away by escaping neutrinos. This 
causes the shock to stall and become an accretion shock before it can drive 
off the envelope of the star |2][3]l^ll5ll6] . 



In the present supernova paradigm, the intense neutrino flux emerging from 
the proto-neutron star heats the matter just behind the stalled shock, eventu- 
ally reenergizing it sufficiently to drive off the envelope and produce a super- 
nova explosion [7f8] . Unfortunately, despite significant progress in supernovae 
modeling, complete self-consistent simulations often fail to produce explosions 
PI5|^PI]pT| . (However, recent results have been more promising [T^13fl4|[T5] .) 
One potential improvement to these models is replacement of incomplete or 
inaccurate treatments of the wide variety of nuclear and weak interaction 
physics that are important in the supernova mechanism. In particular, elec- 
tron captures play a dominant role during the collapse as they significantly 
alter the lepton fraction and entropy of the inner core. These quantities, in 
turn, determine the structure of the core, and the strength and location of the 
initial supernova shock. As a result, the treatment of electron captures signif- 
icantly influences the initial conditions for the entire post-bounce evolution of 
the supernova. 
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Following the pioneering work of Fuller, Fowler and Newman [T6|17fl8|IT9] . 
extended sets of electron capture rates have been calculated for heavy nuclei 
based on rather sophisticated nuclear structure models and advanced com- 
puter algorithms. The rates of nuclei in the mass range of A = 45-65 (around 
100 nuclides; hereafter the LMP pool) are derived from large-scale shell model 
diagonalizations [20|21|22j performed in the complete pf shell at a truncation 
level which guarantees the virtual convergence of the level spectra and the 
Gamow- Teller (GT) strength distributions. In [20|21|22j the rates for electron 
captures (and other weak processes) were derived solely on the basis of al- 
lowed transitions as contributions from forbidden transitions can be neglected 
for those stellar conditions (presupernova evolution) where nuclei in the mass 
range A = 45-65 dominate the matter composition. 

Electron capture rates for about 80 heavier nuclei in the mass range of A = 66- 
112 have been derived within the framework of a hybrid model [23] (hereafter, 
the LMS pool), with the rates calculated via the Random Phase Approxi- 
mation (RPA) based on average thermal nuclear states characterized by oc- 
cupation numbers calculated within the Shell Model Monte-Carlo (SMMC) 
approach pif23] . The SMMC model considers relevant finite-temperature ef- 
fects and correlations among nucleons [22] which both have been identified as 
important for the description of stellar electron capture rates for heavy nuclei 
[23]. Allowed (i.e. GT) and forbidden transitions have been considered in the 
rate calculations of Ref. |24j . 

The consequences of the LMP rates in the presupernova evolution of a massive 
star were studied by Heger et al [27|28] . These rates have also been employed 
to study thermonuclear supernovae (see, e.g., [29|30] ). The effects of improved 
nuclear electron capture during core collapse was studied in Refs. [3T|l32] . 
where the LMSH tabulation was constructed by folding the LMP and LMS 
rates (supplemented by the rates of Fuller, Fowler and Newman for lighter 
nuclei with mass numbers A < 45) with a detailed calculation of the nuclear 
composition assuming Nuclear Statistical Equilibrium (NSE). The assumption 
of NSE among free nucleons and discrete nuclei is well justified during core 
collapse until densities of order 10^'^ g/cm^ are reached. 

The use of the LMP+LMS electron capture rates in supernova simulations for 
a wide range of relevant collapse conditions [3T|32] showed the importance of 
a correct treatment of these nuclear processes. Previous simulations assumed 
that captures on nuclei with neutrons above the p/-shell closure {N > 40) 
were negligible due to Pauli blocking of the GT transitions. (Calculations of 
thermal unblocking of the GT transitions were reported in [53E3], but found 
rather small effects at densities before neutrino trapping. Very recently first 
attempts have been reported to derive the stellar electron capture rates con- 
sistently based on the finite temperature Random Phase Approximation [35] 
or the thermofield dynamics formalism [3H]- Comparison to the SMMC results 
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indicate that the GT unquenching across the shell gap requires higher-order 
correlations than currently accounted for in these approaches [36].) In these 
simulations, electron captures on free protons dominated. However, thermal 
excitation of nuclei and, more importantly, correlations among the nucleons 
can effectively unblock the GT transitions in neutron-rich nuclei, resulting in 
rates which are several orders of magnitude larger than assumed before |26] . 
In addition, although the captures on protons remain faster per particle, the 
total reaction rate for heavy nuclei is greater than for free protons under the 
collapse conditions PT|5^ due to the large abundances of heavy nuclei. As a 
further result, the spectra of neutrinos emitted by the electron capture pro- 
cess are significantly changed as the heavy neutron-rich nuclei have noticeably 
larger Q values than protons. 

An important consequence of the dominance of electron captures on nuclei 
over free protons is the faster decrease of the lepton fraction at high densities 
(p ~ 10^^ g/cm'^) and temperatures (T ^ 1 MeV) as under these conditions the 
nuclear composition is dominated by nuclei for which electron captures were 
originally neglected. Once neutrino trapping sets in at densities around p ~ 
lO'^^ g/cm^ the deleptonization is hindered by final-state neutrino blocking. 
Recently inelastic neutrino-nucleus neutral current reactions have been found 
to have very small effects on the supernova dynamics, but a larger effect on 
the high energy tail of the neutrino spectrum [37] . 

In spite of a smaller inner lepton fraction, resulting in a smaller initial proto- 
neutron star and a weaker shock, the shock ultimately travels slightly further 
out before it stalls in simulations which consider electron captures on nuclei 
[32] . The captures on nuclei from the LMP pool dominate in the outer regions 
of the core. Since these rates are in general smaller than the previously used 
rates from Fuller, Fowler and Newman [TB|17|18|ITU] . the deleptonization in 
the outer regions is smaller. The larger electron fraction in the outer layers 
of the core increases the degeneracy pressure and slows the collapse. A slower 
collapse reduces the rate at which density increases and, hence, the ram pres- 
sure that opposes the shock. In spherically symmetric simulations, this more 
than counteracts the effects of the weaker, deeper initial supernova shock. 

The calculation of supernova-relevant electron capture rates requires an ap- 
propriate nuclear model to determine the individual capture rates as well as 
a reliable account of the many nuclei present in the matter composition and 
of their individual abundances. A shortcoming of the LMSH tabulation was 
the decreasing fraction of the nuclear composition included in the LMP-fLMS 
set of electron capture rates as collapse continues. In fact at the densities 
around neutrino trapping the nuclear composition becomes dominated by nu- 
clei heavier than those considered in this set. To minimize this problem, the 
LMSH tabulation derives an average electron capture rate per heavy nucleus 
on the basis of those nuclei for which individual rates are available. The total 
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rate is calculated from this average rate and the heavy nucleus abundance 
provided by the equation of state. However, this average becomes increasingly 
dominated by a few nuclei, normally the heavier, making it increasingly un- 
certain. It was suggested [21] to use a parameterized electron capture rate for 
the missing nuclei at higher densities. However, a recent study revealed that 
the parameterization is too simple and cannot be used at all conditions during 
the collapse [38] . 

This paper presents electron capture rates for a pool of nuclei which has been 
enlarged in two ways. The first addition is achieved by calculating the capture 
rates within the hybrid model as used in [31] for 170 additional nuclei. The 
SMMC+RPA pool (SMMC pool, below) now covers about 250 nuclei in the 
mass range of A = 66-120 and the proton number range oi Z = 28-45. As the 
SMMC calculations are rather time consuming, it is prohibitive to perform the 
required studies for the many nuclei which are present in the latest stages of 
the collapse. This prompts a new approach to expand the pool further, using 
the fractional occupation numbers of the various shells in the ground state of 
the parent nucleus as calculated from a Fermi-Dirac (FD) parameterization 
in place of the SMMC results. The reduced cost of this approach allows the 
inclusion of more than 2200 additional nuclei to the pool. To improve this 
approach, the parameters of the Fermi-Dirac distribution were adjusted so 
that the electron-capture cross sections matched those available for the SMMC 
pool of nuclei. The Q values needed to calculate the cross sections were derived 
from the predicted masses of the finite-range droplet model jHH], if they were 
not available from Audi et al compilation jlD]. The nuclei added in this way 
are in the range of proton and neutron numbers Z = 28-70 and = 40-160. 
We will refer to them as the FD-I-RPA pool. We further include the FFN 
rates [T6|17|18|IT9] . for nuclei with mass numbers A < 45. Since the majority 
of these 100 nuclei are from the sd shell, we refer to them as the sd pool. 

To have an overview of the pools. Fig. [T] shows the nuclei included in our rate 
evaluation distinguishing the method used for the individual rates. With four 
approaches all contributing, validation is an important consideration. Ideally, 
one would validate the various approaches against experimental data. Unfor- 
tunately, such a procedure is severely limited by the lack of data for excited 
states, while the relative weight of the ground state in the thermal ensemble 
at temperatures present in the collapsing core is rapidly decreasing. Further- 
more, the important role played by the phase space has to be recognized as 
well. The electron chemical potential grows significantly faster than the aver- 
age Q- values of nuclei present in the core. As a consequence, the stellar electron 
capture rates are sensitive to detailed Gamow- Teller distributions only at low 
densities. At higher densities it suffices that the total GT strength and its en- 
ergy centroid are well described. This calls for more elaborate nuclear models 
to be used to derive the capture rates for nuclei at low densities than is re- 
quired for the rates needed at larger densities. We have considered these facts 
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Fig. 1. (Color online) Nuclei included in the calculation of the NSE-averaged rates 
and spectra. The sd pool is marked by circles, the shell model pool is marked by 
plusses, the SMMC+RPA pool is marked by crosses, and the FD+RPA pool is 
marked by diamonds. 



in our strategy for the validation process. At low densities, which is also at 
low temperatures when nuclei are present, the ground states contribute rela- 
tively strongly to the capture rates. For such conditions the core composition 
is dominated by nuclei with the mass numbers A = 45-65 for which the shell 
model diagonalization calculation can be performed. Such shell model rates 
are adopted here. It has been proven by extensive comparison to experiment 
that the shell model describes measured GT+ strength distributions in this 
mass range very well and also gives a very good account of the spectra at low 
excitation energies [20|21|I22|I^ . By the time nuclei with A > 65, for which 
shell model diagonalization calculation are prohibitive due to computational 
restrictions, dominate the core composition, the density, and accordingly the 
electron chemical potential, has grown sufficiently that the capture rates are 
mainly sensitive to the total GT strength and its centroid. For such nuclei we 
have performed RPA calculations with the occupation numbers determined 
from SMMC calculations or from the simple FD parametrization. This pro- 
cedure is validated by comparing the shell model and SMMC+RPA results 
for selected nuclei, indeed finding good agreement, if the electron chemical 
potential is appropriately large. Finally we derive a simple parametrization 
of the SMMC occupation numbers and show that this FD parametrization 
reproduces the SMMC+RPA capture rates quite well. Such a comparison of 
the rates in the overlapping regions together with the outlined reasoning on 
the rate sensitivity to the conditions gives us confidence that the approach 
used is valid. 
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Having produced the individual rates, we have derived the NSE-averaged elec- 
tron capture rates for the appropriate conditions during the collapse based on 
the pool of nearly 2700 nuclei. An extended tabulation of the capture rates and 
the corresponding neutrino spectra is available in electronic form upon request 
from the authors for a wide grid of stellar conditions (defined by temperature, 
the matter density and the electron-to-baryon ratio). This tabulation is most 
appropriate for the study of core collapse super novae, where the large range of 
density and the neutrino spectral information are necessary. This tabulation 
is also appropriate for use in thermonuclear supernovae, where the range of 
density and electron fraction are a subset of those occurring in core collapse 
supernovae, however the value of the wider nuclear range and neutrino spectral 
data considered here is lost. Thus for the thermonuclear supernova problem, 
this tabulation is essentially equivalent to that of Seitenzahl et al. [12] , which 
also folds LMP reactions over an NSE abundance distribution. This tabulation 
is wholly inappropriate for the study of electron capture in X-ray burst ashes 
(see, e.g. [13]), because the ash temperature is insufficient to justify the use of 
NSE. 

The paper is organized as follows. Section [2] describes the hybrid model and 
some computational details. Screening corrections to the rates are also intro- 
duced there, but the corresponding formalism is presented in an appendix. 
Section [3] discusses the obtained results. We conclude with section HJ 



2 Theoretical model 

The hybrid SMMC-I-RPA model was proposed in [2B] to compute electron 
capture rates on nuclei which required such large model spaces for which 
diagonalization shell model calculations are not yet feasible. The hybrid model 
is computationally feasible, but simultaneously incorporates relevant nuclear 
structure physics as configuration mixing (caused by nucleon correlations) 
and thermal effects. A pairing-|-quadrupole residual interaction [33] was used 
which avoids the sign problem associated with using realistic interactions in 
SMMC studies [2S]- SMMC calculations at finite temperature are used to 
obtain occupation numbers for the various neutron and proton valence shells 
in the parent nucleus, which are then used to calculate electron capture cross 
sections and rates within a Random Phase Approximation (RPA) approach 
with partial shell occupancies (the method is explained in Ref. [3S])- 

We have performed hybrid SMMC-I-RPA calculations of electron capture rates 
for additional nuclei extending the pool of nuclei used in [23|3T] to more 
neutron-rich nuclei and filling occasional holes in isotopic chains. As mentioned 
above, the first step in the hybrid approach is to obtain the shell occupation 
numbers for protons and neutrons using the SMMC. Using the same interac- 
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tion and model space as in [31] (i.e. full pf-gds shells with 50 valence orbitals 
for protons and neutrons), we have calculated such occupation numbers at 
finite temperatures for 168 nuclei with < 61. For nuclei with even larger 
neutron number we switched to the f5/2P-gds-hii/2 model space, changing the 
pairing and quadrupole strength parameters of the residual interaction to ac- 
count for the change in the model space. These parameters were adjusted such 
as both model spaces predict similar properties for selected nuclei. The single- 
particle energies were taken from |16]. The energy of the /iii/2 was determined 
using the same Woods-Saxon parameters as in placing it at 15.12 MeV 
above the f'i/2 shell. 

The radial wave functions as well as single-particle energies for the RPA cal- 
culation (i.e. the second step for the hybrid method) were taken from the 
Woods-Saxon potential: 
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The depth of the potential was adjusted to reproduce the proton separation en- 
ergy in the parent nucleus and the neutron separation energy in the daughter 
nucleus, as discussed in [31]. Other parameters of the Woods-Saxon poten- 
tial were: tq = 1.27 fm, R = tqA^^^, a = 0.65 fm, A = 32. The separation 
energies were calculated using masses from the Audi et al compilation |40j . 
supplemented by predictions from the finite-range droplet model [32] • The 
same masses were used to evaluate the Q value for the calculation of electron 
capture cross sections as will be discussed later. 



Following the spirit of RPA+BCS calculations we have shifted the energies of 
the unoccupied proton states by the amount of the pairing energy estimated 
as 12/\/A MeV to avoid the appearance of spurious zero energy transitions 
in situations where orbits are partially filled. It should be noted that this 
energy shift does not affect the physical transitions that take place between 
occupied proton states and unoccupied neutron states. Ref. [31] has used a 
constant energy shift of 2.5 MeV which leads to slightly larger capture rates. 
For consistency we have repeated the calculations of the capture rates for the 
LMS pool of nuclei with the A-dependent energy shift. 

In the RPA calculations we used the Landau-Migdal force as the residual inter- 
action with the parameters taken from [17] , except for the overall scale factor 
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Fig. 2. (Color online) A comparison of the electron capture rates on ^^'^^Ni cal- 
culated from the diagonalization shell model (only allowed contributions) and the 
hybrid SMMC+RPA model (both allowed and forbidden contributions). Stellar con- 
ditions of the 25Mq trajectory (see Tabled]) are used. 

Cq = 302 MeV fm^ [JH]. In the calculation of the rates we included all multi- 
pole transitions with J < 3. The dependence of the operators on momentum 
transfer has been considered which leads to a reduction of the capture cross 
sections at high electron energies [23]. Motivated by global RPA calculations 
for muon capture on nuclei [19|50] . the J = 1^ multipole transitions have been 
quenched by a factor of (0.7)^ but not the other. 



The diagonalization shell model has been clearly established as the method of 
choice to calculate electron capture rates, provided the required model space 
is feasible with presently available computer memory. The later restriction, 
however, prevents its use for many of the heavier neutron-rich nuclei that 
are present in the collapsing stellar core, forcing us to use the hybrid model 
as discussed above. We have tested the validity and the accuracy of this hy- 
brid SMMC-I-RPA approach by comparing electron-capture rates, calculated 
within this model, against diagonalization shell model results for two nickel 
isotopes (Fig. [2]). The calculations have been performed for a set of stellar 
conditions (density, temperature. Ye value) appropriate for the inner core of a 
collapsing 25Mq star (see Tabled]). The shell model rates are based only on 
the allowed Gamow- Teller (GT) transitions, while the SMMC-I-RPA calcula- 
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Fig. 3. (Color online) A comparison of the neutrino spectra emitted during electron 
capture on ^^'^^Ni calculated from the diagonalization shell model (only allowed 
contributions) and the hybrid SMMC+RPA model (both allowed and forbidden 
contributions). Stellar conditions of two zones (10 and 15) of the 25Mq trajectory 
(see Table [I]) are used. 

tion also includes forbidden transitions. To identify the effect of such forbidden 
transitions we have also performed SMMC+RPA calculations considering only 
GT transitions. As can be seen in Fig. O the shell model and the hybrid rates 
are quite similar for electron chemical potentials fie larger than about 10 MeV 
when only the Gamow- Teller transitions contribute to the electron-capture 
rate. At lower fi^ values the capture depends more sensitively on the detailed 
structure of the GT strength distribution which is better described by the 
diagonalization shell model than in the hybrid model. As explained in [5T] . 
with increasing electron chemical potential the capture rate becomes more and 
more dependent solely on the total GT strength rather than on the details of 
its distribution. It is fortunate that conditions that produce such low electron 
chemical potentials, for which a detailed reproduction of the GT strength 
is required for a reliable description of the rate, also result in relatively low 
mass nuclei in NSE, for which diagonalization shell model calculations can 
be performed and which are included in the LMP pool. Forbidden transitions 
start to contribute noticeably to the capture rates for fi^ > 30 MeV. At such 
conditions corresponding to rather high densities the nuclear composition is 
dominated by nuclei heavier than included in the LMP pool and hence the 
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Fig. 4. (Color online) Single j'-shell fractional occupation numbers in 32^6 at dif- 
ferent SMMC temperatures. Proton values are shown on the left, neutron val- 
ues are shown on the right. The fractional j-shell occupancy is calculated by 
Uj = nj/{2j + 1), where rij is the number of particles in the shell. The single- 
particle energies are taken from the nucleus- adjusted Woods-Saxon potential. 
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Fig. 5. (Color online) Single j-shell fractional occupation numbers for 28^1 and 
^SjyRb as a function of single-particle energies. The SMMC results are shown by 
diamonds; proton values are marked by full symbols, neutron values are marked 
by empty symbols. The values for 28-'^i and ^-jyRb were obtained at temperatures 
T = 0.50 and 1.33 MeV, respectively. The fractional occupations given by the 
Fermi-Dirac distribution for four values of Tpo are shown by lines: Tpo = 2.50 MeV 
(dotted), TpD = 1.67 MeV (dash-double-dotted), Tpo = 1-25 MeV (dash-dotted), 
Tfd = 1.00 MeV (dashed). 

neglect of higher multipole transitions in the LMP rates is unimportant. 



The reasonable reproduction of the shell-model Gamow- Teller distribution by 
the SMMC+RPA approach is also reflected in the emitted neutrino spectra, 
shown in Fig. [3l Once the value of the electron chemical potential gets suf- 
ficiently large, the emitted neutrino spectra obtained within the shell model 
and SMMC-I-RPA approaches are quite similar. 

A comment about the occupation numbers used in the SMMC-I-RPA cal- 
culations is in order. In principle one should calculate these numbers at all 
the temperatures needed to construct the stellar capture rate table. However, 
such a procedure is computationally untenable. Fortunately it turns out that. 
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within the range of temperatures relevant to the phase of the collapse at which 
electron captures are important (T ^ 0.8-1.4 MeV), the occupation numbers 
do not vary too much. This is demonstrated in Fig. H] which shows the proton 
and neutron occupation numbers for ''^Ge obtained by performing SMMC cal- 
culations at different temperatures. We conclude from this comparison that 
it might suffice to calculate the SMMC occupation numbers only at a single 
temperature which we chose as the value where the respective nucleus has a 
large relative weight in the NSE composition. 

SMMC calculations are rather time-consuming. Thus calculations to extend 
the SMMC pool beyond the 250 nuclei considered here to include the couple 
of thousand species which are present in the supernova composition is pro- 
hibitive. However, the observation that the SMMC occupation numbers do 
not vary too much with temperature motivated us to find a parameterized 
form for the SMMC occupation numbers which could then readily be used 
in RPA calculations of electron capture rates. This goal is achieved by as- 
suming a Fermi-Dirac parameterization for the proton and neutron fractional 
occupation numbers of various shells with energy e: 

i^fd(£^,7fd) = — 77 rjTf — r- (2) 

1 + exp{{e - ;u)/Tfd| 

The chemical potentials are fixed by the total proton and neutron numbers 
(here denoted by A^) 

N = J2{2j, + l)uF^{e,TFB). (3) 



Here jg denotes the total angular momentum of a shell having the single- 
particle energy e. These energies are taken from a Woods-Saxon potential, the 
depth of which is adjusted such as the neutron and proton chemical potentials 
resulting from equation equal the respective neutron and proton separation 
energies, taken from either experiment or the compilation of ref. [32] • The 
number of the considered single-j shells depends on the numbers of nucleons 
and the value of the temperature Tpo and exceeds that used in the SMMC 
studies. The Fermi-Dirac distribution has one more undetermined parameter, 
TpD- We will vary it as a parameter below and fix its value by attempting to 
reproduce the SMMC+RPA electron capture cross sections for the 250 nuclei 
present in the LMS pool. In general we have included at least 2 major oscillator 
shells in our RPA calculations outside a closed N = Z = 20 core. 

Fig. [5] compares the SMMC occupation numbers for two different nuclei with 
those obtained by an FD parameterization for different values of the parameter 
TpD- The values Tpo = 1-0 MeV and 1.25 MeV give a rather fair reproduction 
of the SMMC occupation numbers for the two nuclei shown. (As we will see 
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in the next chapter the SMMC+RPA capture cross-sections are globally best 
reproduced by the value Tfd = 1-67 MeV.) It is important to note that the 
occupation numbers reflect both correlation and thermal effects. Thus, unlike 
for a system of non-interacting particles, the parameter TpD should not be 
interpreted as the temperature of the nucleus. 

The stellar electron capture rate Aec on a particular nucleus is related to the 
electron capture cross-section dec by: 

1 °° 

■^^^ = J Pl'^ocise, Si, ef)f{ee, /ie, T)dee (4) 

if -0 



where 6° = max(Qj/-, mgC^), Pe = {si — mlc^Y^"^ / c is the momentum of the 
incoming electron with energy Se and me the electron rest mass. Under the 
conditions present in the collapsing core of a supernova, electrons obey a 
Fermi-Dirac distribution /(ee,/ie,T) = [1 + expK^g — l^'e) / k^T}]'^ with tem- 
perature T and electron chemical potential /^e- 0"ec(^e5 Si,Sf) is the cross section 
for capture of an electron with energy from an initial proton single particle 
state with energy Si to a neutron single particle state with energy Sf. The 
cross section is computed within the Random Phase Approximation as de- 
scribed above. Due to energy conservation, the electron, proton and neutron 
energies are related to the neutrino energy, e^, and the Q- value for the capture 
reaction [51] : 

Qif — Se Sy = £ ^ 5^ , (5) 

e--e^^=e*. + fi + /\np. (6) 



Here we take jl = fin — fip, the difference between neutron and proton chemical 
potentials in the nucleus and A„p = (M„ — Mp)c^ = 1.293 MeV, the neutron- 
proton mass difference. 

Equation (jSD constitutes the definition of the quantity e*j. At zero temper- 
ature e*j corresponds to the excitation energy in the daughter nucleus. For 
the transition from the initial ground state to the daughter ground state the 
excitation energy must be zero. This fact is used for fixing the value of fi: 

Qoo = {Mf - M^y = fi + A„p, (7) 



where Mj and Mf are the nuclear masses of the parent and daughter nuclei 
in the electron capture reaction, and the values of /ip and selected to 

equal proton and neutron separation energies, S'p(Z, A^) and S'„(Z — 1, A^+ 1). 
At finite temperature the value of e*j can be negative when the nucleus is 
de-excited by capturing an electron. 
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The capture rate for individual nuclei is affected by screening corrections that 
are implemented in our calculations as discussed in the appendix. The screen- 
ing corrections affect the electron capture in two ways: they effectively increase 
the Q value and lower the electron chemical potential. As discussed in the next 
section and the appendix, both effects lead to a reduction of the electron cap- 
ture rate in the medium as compared to the undisturbed case. 

Calculation of the partial cross sections allows us to obtain the spectrum of 
the emitted neutrinos. To have the spectrum normalized to unity, we define it 
as follows: 

1 1 

ec f 



with Be = e f — Si + Bp and pe the corresponding electron momentum. Equa- 
tion dH]) defines the spectrum of emitted neutrinos as produced by captures on 
a particular nucleus. 

During stellar collapse, the matter is composed of individual nuclei until den- 
sities of order 10^'^ g/cw? are reached. While presupernova studies (i.e. simu- 
lations which cover the late-stage evolution of a massive star until the inner 
core has reached densities up to 10^ g/cm^) require explicit consideration of 
the extensive nuclear networks, temperatures in the subsequent supernova col- 
lapse phase are high enough to bring reactions mediated by the strong and 
electromagnetic force into equilibrium with their inverse reactions and hence 
the nuclear composition can be described by Nuclear Statistical Equilibrium 
(NSE). It is customary to split the required electron capture rates for super- 
nova simulations into rates for protons Rp = YpXp, and for nuclei Rh = J2i '^i^i 
where the index i runs over all isotopes other than protons present in the stel- 
lar composition. The required abundances for protons Yp and for the various 
nuclear species Yi are derived from a Saha-like NSE distribution of the indi- 
vidual isotopes where we included consideration of degenerate nucleons and 
plasma corrections to the nuclear binding energy using the screening for- 
mula of Slattery et al [53] with the parameters from [51] and the extension 
to the weak screening regime as suggested by Yakovlev and Shalybkov [55] as 
discussed in the appendix. Compared to the LMSH tabulation used in [3T|32] , 
we have extended the pool of nuclei in the NSE distribution to heavier nuclei 
(Ag to At). Furthermore we have used recently published partition functions 
which extend to temperatures in excess of 10 GK [56]. The electron chem- 
ical potential necessary for the calculation of the different electron capture 
rates are determined using the helmholtz Equation of State by Timmes and 
Swesty [57] . 

In [HT|^ the stellar electron capture rate (Aec) has been derived as the NSE- 
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average over all nuclei for which individual rates are available. Hence, 

(Aec) = (9) 



Similarly, the NSE-averaged neutrino spectrum is obtained by calculating the 
ratio 

= JKT) = E.5^Af 



This spectrum is normalized to unity, with the absolute neutrino emission rate 
being (Aec) (AC (E,)). 

However, there is a potential problem with the average rate as defined above, 
when applied to the pool of nuclei used in [3T|32] . as this pool was limited to 
nuclei with mass numbers A < 112, sampling only a rather small fraction of 
the total NSE abundance in the later phase of the collapse (changes in nuclear 
composition during the collapse are illustrated in [58] )• While this averaging of 
the rate prevents the total rate of electron capture from unphysically dropping 
to zero as the fraction of nuclei with calculated rates declines, the average rate 
becomes increasingly dominated by the species at the neutron-rich edge of the 
calculated pools, making the averaged rates of pill^ uncertain. It is the aim 
of the present work to overcome this shortcoming by appropriately enlarging 
the pool of nuclei from which the averages iQ and ffTOj) are being derived. The 
present pool contains nearly 2700 nuclei (in contrast to about 280 considered 
in [3T|32] ). For nuclei with mass numbers A < 65 our pool consists of the 
LMP shell model rates [22] (i.e. the p/-shell nuclei) supplemented by the FFN 
rates for the lighter sd-shell nuclei. We have also extended the pool of nuclei 
for which rates are based on the hybrid SMMC+RPA approach to include 
nuclei up to A = 120 (SMMC pool). The pool is completed by more than 
2500 nuclei for which we have derived rates based on RPA calculations with 
the occupation numbers approximated by the parameterized FD distribution 
as discussed above (FD pool). As some of the nuclei calculated within this 
approach are overlapping with the nuclei from other pools, we do not use the 
FD+RPA results if the rates and spectra are already provided by other pools, 
unless specifically indicated. 

In the next section we present electron capture rates and the emitted neutrino 
spectra for conditions during the core collapse of two progenitor stars with 
different masses. The first set of conditions follows a central mass element 
(at an enclosed mass of 0.05 solar masses) during the collapse of a 15M0 
progenitor star derived from a simulation using electron-capture rates from 
the LMSH tabulation [3T|32] and general relativity. The second set is based 
on a similar simulation, but for a 25M0 progenitor star. The stellar conditions 
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Table 1. Conditions for a central mass element in two progenitor stars at several points in time during the collapse. A and Z is the 
average mass and proton numbers as given by the NSE. Yp is the abundance of protons. Yf^ is the abundance of heavy nuclei defined as 
Yh = Yi) — Yp — Yn — Ya, i.e. everything except for protons, neutrons and alpha particles. Zh and Ah are the average mass and charge 
of the heavy nuclei. The last two columns list values of the unscreened and screened electron capture rates ((A°^) and (Aec), respectively). 
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Fig. 6. (Color online) Fraction of nuclei covered by the various pools of nuclei as 
defined in the text. The fractions have been calculated for the two stellar trajectories 
given in Table [D The pools are sd (dotted hne), LMP (dashed hue), SMMC+RPA 
(double-dash-dotted hne), and FD+RPA (dash-double-dotted line). Sohd hues show 
the summed pool coverage. Thick solid lines show present pool coverage, and thin 
solid lines show coverage by the LMSH pool. (^^ ^i)nucici is calculated by summing 
over all nuclei except protons, neutrons and a particles. 

are defined in Table [H 



3 Results and discussion 

As mentioned above, the pool of nuclei adopted to derive the LMSH elec- 
tron capture tabulation used in [3T|32] was limited to isotopes with A < 112 
and turns out to cover only a small fraction of the total NSE abundance in 
the latest stages of the collapse for two reasons. First, heavier nuclei appear 
in the nuclear composition with noticeable abundances. Secondly, the inclu- 
sion of electron capture on nuclei makes the matter more neutron-rich than 
anticipated in prior studies which considered only capture on protons and ne- 
glected the dominant electron reducing weak process. As the pool of nuclei 
used in [31] has been constructed based on trajectories from prior collapse 
simulations which solely considered capture on protons, the pool was missing 
relevant neutron-rich isotopes. 

The limitation of the nuclear pool used in (LMSH pool, i.e. the combined 
sd, LMP and LMS pools) is shown in Fig. El This sample is complete for the 
early stage of the collapse (< 10^*^ g/cm^) where the pool includes more than 
90% of the nuclei present in the supernova medium (indices 1-3 of the ISMq 
and 1-4 of the 25Mq trajectories in Table [1]). However, the coverage drops 
quickly as the collapse progresses. The pool covers more than 50% of the total 
NSE abundance for densities < 10^^ g/cm^ (indices 7 and 10 respectively). At 
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Fig. 7. (Color online) Electron capture rates for four nuclei from the SMMC pool: 
75,76 2n, iii.ii2]Nji3^ r^^ie SMMC+RPA rates are shown by sohd lines. The FD+RPA 
rates have been calculated for TpD = 2.50 MeV (dotted line), T^d = 1-67 MeV 
(dash-double-dotted line), Tfd = 1-25 MeV (dash-dotted line), and Tfd = 1-00 MeV 
(dashed line) . The stellar conditions are those of 25Mq progenitor star (Table [T]) . 

the conditions around neutrino trapping (indices 10 and 12, respectively, of the 
15Mq and 25Mq collapse trajectories) only 20% of the total NSE abundance 
is sampled by the LMSH pool. This drastic underrepresentation is improved 
slightly by including the additional 170 nuclei for which we have calculated 
capture rates using the SMMC+RPA approach, pushing the limits at which 
the 20% coverage of the total NSE abundance is reached to indices 11 and 15 
for the 15Mq and 25Mq trajectories, respectively. 



Obviously improving the coverage requires an appropriate enlargement of the 
nuclear pool from which the capture rates are derived. This is in particular 
motivated by the fact that the omission of the neutron-rich nuclei from the 
pool in [H] could result in a systematic overestimate of the rates as the capture 
process gets increasingly hindered with growing neutron excess due to Pauli 
blocking. However, covering the relevant nuclei by SMMC+RPA calculations 
is too computationally demanding. Thus we have performed evaluations of the 
capture rates for more than 2200 nuclei using occupation numbers derived from 
a FD parameterization. These nuclei consider the isotope chains for charge 
numbers Z = 28-70. When these nuclei are added to the pool. Fig. [B] shows 
that at least 60% of the total NSE abundance is considered during the entire 
collapse evolution until densities of order 10^^ g/cm^ are reached and the 
description of the nuclear composition by individual nuclei is inadequate. We 
stress that at the beginning of the core collapse, the center of a progenitor star 
is populated by iron group nuclei, therefore electrons are mostly captured by 
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Fig. 8. (Color online) NSE- averaged electron-capture cross-sections as calculated by 
the SMMC+RPA and FD-hRPA approaches for the 250 nuclei of the SMMC+RPA 
pool at three stellar conditions of the 25Mq trajectory: the snapshot number 4 is 
used in the left panel, the number 10 in the center panel, and the number 13 in the 
right panel. The notation is the same as in Fig. [71 
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Fig. 9. (Color online) NSE-averaged electron-capture rates as calculated by the 
SMMC-FRPA and FD-^RPA approaches for the 250 nuclei of the SMMC-^RPA pool. 
The stellar conditions are those of 25Mq progenitor star (Table [1]). The solid line 
shows the SMMC-I-RPA rate. The dashed line gives the FD-I-RPA rates calculated 
with Tfd = 1.67 MeV. 




the nuclei from the LMP pool. As the collapse progresses, the nuclei become 
heavier and more neutron rich. Nuclei from the SMMC+RPA pool dominate 
the rates at densities from 5 x 10^° g/cm^ to neutrino trapping around 5 x 
10^1 g/cm^. At even higher densities nuclei from the FD+RPA pool dominate 
the rates. 



The FD+RPA approach is considered as an approximation to SMMC+RPA 
calculations. To verify its validity we previously compared the FD-predicted 
occupation numbers to those obtained from the SMMC calculation and found 
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Fig. 10. (Color online) Pool-averaged electron-capture rates calculated along the 
stellar trajectories for the ISMq and 25Mq progenitor stars. The rates based on 
the sum of all pools of nuclei are shown by solid lines. The dashed lines show the 
average rate when the FD-I-RPA pool is omitted. The dotted lines show the average 
rate for the sum of all pools when the screening effects to the rates are neglected. 

fair agreement (Fig. |5]). The SMMC occupation numbers show some deviation 
from the smooth FD behavior, caused by correlations introduced by the resid- 
ual interaction, and not reproduced by the FD distribution. These deviations 
translate into differences of the rates for individual nuclei (see Fig. [7]). How- 
ever, they do not represent a systematic effect and are smeared out when the 
averaging over many nuclei takes place. This is demonstrated in Fig. [8] where 
we compare electron capture cross sections as derived in the SMMC+RPA ap- 
proach with those obtained in the FD-I-RPA approximation. The comparison 
is made for the 250 nuclei, for which SMMC-I-RPA rates are now available, 
adopting for the NSE averaging the stellar conditions of the 25Mq progenitor 
star (Tabled]). We note that the SMMC pool only includes heavy nuclei with 
charge numbers Z > 28, while the capture on the lighter pf shell nuclei has 
been evaluated on the basis of the diagonalization shell model and is here 
summarized within the LMP pool. Fig. [H] demonstrates the results for the 
cross section comparison at three distinct points along the stellar trajectory. 
The left plot corresponds to a rather early stage of the collapse (presupernova 
phase, index 4) where the SMMC pool covers 81% of the total abundance 
of heavy nuclei with Z > 28. The center plot corresponds to collapse con- 
ditions just before neutrino trapping sets in (index 10). Here the 250 nuclei 
of the SMMC pool contribute 99% of the abundance of heavy nuclei with 
Z > 28 and dominate the capture rate in the combined pool. The right plot 
represents conditions after neutrino trapping (index 13). The coverage of the 
heavy nuclei abundance by the SMMC pool has now decreased to 60% and the 
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Fig. 11. (Color online) Pool-averaged emitted neutrino spectra for the 15Mq and 
25Mq trajectories. The line legend is the same as in Fig. [TOj Two stellar conditions 
are used in each case corresponding to snapshot numbers 10 and 15 of the respective 
trajectory. For snapshot number 10 in the lower panel the curves "no FD" and "full" 
coincide. 

NSE-averaged cross-sections are dominated by the most neutron-rich nuclei 
included in the pool. The FD+RPA cross sections have been calculated for four 
different values of the parameter Tpo- The cross sections at fixed electron en- 
ergy grow with increasing TpD as more nucleons are excited across the = 40 
shell gap opening neutron holes in the pf shell and allowing for proton excita- 
tions within the gds shell. Differences between the calculations with different 
values of TpD are most noticeable at low electron energies where the capture 
process is sensitive to the details of the GT strength distribution. In turn, the 
agreement between the calculations improves for growing electron energy. We 
find that the SMMC-I-RPA results are best reproduced by the NSE-averaged 
FD-I-RPA cross sections for Tpo = 1-67 MeV. We will use in the following 
this value for Tpo- Fig. IH] shows that the FD+RPA approach with TpD = 1-67 
MeV indeed describes the NSE-averaged SMMC electron capture rates quite 
well for chemical potentials fig > 10 MeV. This is the regime of interest as we 
will use the FD+RPA approach to estimate the capture rates for neutron-rich 
nuclei which only contribute significantly to the total NSE abundance at high 
densities, i.e. at the conditions for which the electron chemical potential is 
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larger than 30 MeV. 



The discussion above took advantage of the extended SMMC+RPA pool of 
250 nuclei, which is noticeably larger than the one adopted in pT|32] . This 
increase in coverage lowered the NSE-averaged rates [38] at higher values of 
electron chemical potential as many of the added nuclei are more neutron- 
rich than those used in the original LMSH rates [31]. As explained above, an 
increase in the neutron excess usually leads to larger Pauli blocking of GT 
transitions and hence smaller rates. Furthermore, more neutron-rich nuclei 
have larger Q values which also reduces the electron capture rate. 

Finally we compare electron capture rates for different pools of nuclei in 
Fig. [i The srf+LMP+SMMC+FD pool ("all pools" in the figure legend) 
considers around 2700 nuclei, including more than 2200 nuclei for which the 
rates have been derived using the FD+RPA approach. The sd-l-LMP-l-SMMC 
pool ("no FD" in the legend) leaves out the FD+RPA rates. A comparison 
between the NSE averaged rates for these pools allows an estimate of the rel- 
evance of the omission of the heavy and most neutron-rich nuclei for stellar 
electron capture. The rates have been calculated for the conditions along the 
collapse trajectories of the ISMq and 25M0 progenitor stars as listed in Ta- 
ble [U (The obtained rates are also given in this table.) As can be observed 
in Fig. [10] the additional inclusion of the 2200 heavy and neutron-rich nuclei 
lowers the NSE-averaged rates at larger electron chemical potential slightly, 
where the added nuclei dominate. 

Fig. [To] also quantitatively demonstrates the impact of the medium effects on 
the NSE-averaged rate. As mentioned above, screening effects on the electron 
capture rates lead to a reduction of the electron capture rates which can 
amount to almost a factor of 2 at large densities (large chemical potentials 
in Fig. [TU]) . This is in addition to the effects of screening on the strong and 
electromagnetic rates that determine the NSE [52][5^ . which is included in 
[3T]|32] . Despite this modest reduction in the rates electron capture on nuclei 
still dominates over capture on free protons during the stellar collapse. 

As mentioned in the first section, a very important quantity for supernova 
simulations is the spectral distribution of the emitted neutrinos which becomes 
relevant especially at higher densities when neutrinos are becoming trapped. 
Fig. [11] illustrates the calculated spectra of the neutrinos emitted at different 
stages of the collapse using the conditions defined by indices 10 and 15 of the 
collapse trajectories of the 15M0 and 25Mq progenitor stars (Table [T]). The 
neutrino spectra are only slightly changed if we increase the pool of nuclei by 
more than 2200 heavy and neutron-rich nuclei. At the beginning of the collapse 
these nuclei do not contribute (this is especially true for the snapshot number 
10 of the 25Mq star, where the lines "full" and "no FD" are indistinguishable 
in Fig. [lU). At /ie ^ 20 MeV, the SMMC+RPA nuclei dominate. For the later 
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collapse phase the inclusion of the FD+RPA nuclei slightly reduces the average 
neutrino energy, in agreement with the fact that the heavy neutron-rich nuclei 
have slightly larger Q values. 

Already during the construction of the LMSH tabulation [31] it was realized 
that the shell model approaches can provide rates only for a sample of the 
nuclei present in the medium of a collapsing star. For nuclei not included 
in the sample, it was suggested to use a formula based on phase space con- 
siderations with parameters fixed to approximately reproduce the rates of 
nuclei in the pool (eq. (1) in [31]). A more careful study of this proposal re- 
vealed a deficiency in this approach [38] . For electron chemical potential above 
/Xe > 40 MeV, forbidden transitions contribute significantly, and the value of 
the "typical" transition matrix element 5, derived for pure allowed GT transi- 
tions, becomes too small. We do not suggest an alternative parameterization of 
the rates here. However, we would like to exploit the observation that electron 
capture rates are becoming insensitive to the detailed structure of the nuclear 
transitions with increasing chemical potential and are then relatively simple 
functions of Q for a fixed yUg. Using the 250 nuclei of the SMMC pool we have 
derived an (unweighted) average rate for fixed values of {Q,fj,e)- Identifying 
each of the more than 2200 nuclei, for which the rate has been calculated solely 
on the basis of the FD parametrization, by its Q value, we have derived an 
alternative set of capture rates where the rate of each nucleus, as determined 
by the FD+RPA approach, has been replaced by the average rate from the 
SMMC pool at the corresponding values of Q and /ig. The two sets of rates are 
compared in Fig. [121 The agreement between the sets is excellent, providing 
confidence that the limitations of FD+RPA approach are of little impact here. 
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Fig. 12. (Color online) Pool-averaged electron capture rates for the trajectory of 
the 25Mq progenitor star. The rates based on the combined sd, LMP and the 
SMMC+RPA pool is shown by the dotted line. The averaged rate obtained when 
the FD+RPA pool of rates is also included is shown by the solid line. The dashed 
line shows the rates if the SMMC-sampled {Q, A) rates (see text) are used to replace 
the FD+RPA results. 
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4 Conclusions 



Core-collapse supernova simulations require reliable electron-capture (and neu- 
trino emission) rates to track the deleptonization process. In this pursuit we 
present here an enlarged pool of nuclei from which we derive NSE-averaged 
electron capture rates and the corresponding emitted neutrino spectra. This 
extended pool now consists of around 2700 nuclei. The rates of the sd shell 
nuclei are taken from Refs. [T6|17fl8|[T9] . modified by appropriate screening 
corrections. The pf shell nuclei are adopted from the LMP rates [22]. The 
SMMC+RPA set [2^1151] was enlarged to include in total 250 nuclei [M]. Ad- 
ditionally, more than 2200 nuclei are calculated within the FD-I-RPA approach, 
adjusted to reproduce the SMMC-I-RPA rates. In the FD-I-RPA approach the 
parent nucleus occupation numbers at finite temperature are approximated 
by a Fermi-Dirac parameterization taking the parameter value Tfd = 1-67 
MeV. The latter three sets of rates include screening effects directly, while 
for the FFN rates we use an approximate, but rather accurate, prescription. 
This huge number of individual electron-capture rates was used to obtain the 
averaged electron-capture rates at various stellar conditions (p, T, Fg)- The 
resulting table is available upon request from the authors. 
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A Screening corrections to electron capture rates 

Coulomb corrections are known to play an important role in determining the 
thermodynamical properties of a high density plasma [55]. At the conditions 
we are interested in this work the matter can be assumed to be in nuclear 
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statistical equilibrium. This regime has been studied by Bravo and Garci'a- 
Senz [59] and we have generalized their treatment to the study of electron 
capture rates. When Coulomb corrections are included the chemical potential 
of the nuclear species is given by: 

/^i = /ii,o + /^i,c- (A.l) 



Here /ij^o is the chemical potential in the absence of Coulomb effects which 
is generally given by the Boltzmann statistics, and fii^c is the contribution to 
the chemical potential due to the interaction of nucleus i with the electron 
background. 

The core of the star constitutes a multicomponent plasma that we will treat in 
the additive approximation. In this case, all the thermodynamic quantities are 
computed as the sum of the individual quantities for each species. If one further 
assumes that the electron distribution is not affected by the presence of the 
nuclear charges (uniform background approximation), the Coulomb chemical 
potential of species i is given by 



fii,c = kBTfciT,), (A.2) 

where fc is the Coulomb free energy per ion in units of ksT and is the 
ion-coupling parameter, 

75/3 2 

r. = ^Pr, = (A.3) 

where is the electron sphere radius, = (3/(47rne))^/'^, with rif, the electron 
density. 

For the free energy in the regime F > 1, we use the expression 

/c(F) = aF + 46Fi/^ - AcT-^/^ + rflnF + e, (A.4) 

with the values of the parameters a, b, c, d and e taken from [61] : 

a = -0.898004, 6 = 0.96786, c = 0.220703 

(^= -0.86097, e = -2.52692. (A.5) 



For F < 1, we use the expression suggested in ref. 

fc{T) = -^T'/' + ^F^. (A.6) 
V3 7 
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The first term reproduces the Debye-Hiickel hmit for F ^ 1 and the parame- 
ters P and 7 are determined requiring the continuity of the internal energy and 
its derivative at F = 1. These conditions give /3 = 0.295614 and 7 = 1.98848. 

The Coulomb corrections contribute to the electron capture rates in two dif- 
ferent ways. First, due to the fact that the chemical potential depends on 
Z, screening effects will change the threshold energy for the capture by the 
amount |62j : 

AQc = fic{Z-l)-fic{Z), (A.7) 



where Z is the charge number of the capturing nucleus. As fic is negative, 
this correction increases the energy threshold and thus reduces the electron 
capture rate. 

Second, the energy of the captured electron is affected by the presence of 
the background electron gas. Its energy will be reduced compared to the un- 
screened case. The magnitude of this effect can be determined using linear 
response theory [63]. Moreover, we can assume that the screening potential, 
Vs, is constant inside the nucleus with a value which can be determined from 
equation (17) of Itoh et al. [63] evaluated at the nuclear radius. 

To consider both screening corrections the phase space integral [21] gets mod- 
ified. The rate of the electron capture during the nuclear transition from an 
initial state i to a final state / ^ is given by: 

1 °° 

^^f = / Pl^eciSe, ei,€f)f{ee + K, /^e, T)&e, (A.8) 



= max(Q|j, mgC^), and the energy of the emitted neutrino is related to the 
electron energy by = Se — Qfj, where 

Ql^ = + AQc, (A.9) 
with Qif being the capture threshold in the absence of screening corrections. 



Note that here the state denotes either a single particle state if the rates are 
determined using an RPA approach as in the present calculations or a many body 
state if the shell-model or similar approach is used. 
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Assuming a Fermi-Dirac distribution for the electrons and using the property 
/(^e + Vs, /ie, T) = f{ee, /ig — Vg, T) oue can rewrite equation (lA.SP as follows: 



\ec _ 



1 °° 

— J plaec{ee,ei,€f)f{ee,He-Vs,T)dee. (A.IO) 



The presence of the electron background effectively reduces the electron chem- 
ical potential and hence also reduces the capture rate. In summary, in order 
to evaluate the electron capture rates under the presence of screening the 
threshold energy should be modified following equation flA.9p and the chemi- 
cal potential of the electrons should be reduced by Vg. 



The formalism discussed above allows to include the screening corrections in 
the evaluation of electron capture rates provided that one uses some nuclear 
model to determine the initial and final states. However, many currently avail- 
able tabulations of electron capture rates |16II17II18|19|22|64|65|66] have been 
computed without including screening corrections. (Note that for this paper 
we recalculated the LMP rates [22] with the screening corrections included di- 
rectly.) For such tabulations it is thus important to determine an aposteriori 
prescription to incorporate screening effects for the use of such rate tables in 
astrophysical simulations. 

A very good approximation can be obtained by using the effective ft value 
formalism introduced by Fuller et al. [19] which allows accurate interpolations 
of the tabulated electron capture rates. In this approach an effective /t-value is 
defined that represents an average nuclear matrix element which characterizes 
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Fig. A.l. Rate for electron capture on ^^Fe for a temperature of 7.5 GK as a function 
of the electron density, pY^. The solid line corresponds to the rate computed in 
ref. [22], the dashed line uses the same matrix elements but includes screening 
corrections as discussed in the text, see eq. (|A.10p . The dashed-dotted line has been 
obtained using the approximation suggested in equation ()A.13p . 
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the capture rate on a given nucleus: 

/teff = ln2:^^^^, (A.ll) 



where Aec is the tabulated electron capture rate, Q is the ground-state to 
ground-state transition energy, and F is the integral 

1 °° 

/i,) = / w\Q + wffiw, /ie, T)dw, (A.12) 

max(Q,me c^) 



Assuming that the average nuclear matrix element ft^s is unaffected by screen- 
ing, the ratio of electron capture rates with and without screening is equal to 
the ratio of / functions, leading to: 

F{Q + AQc,fie-Vs) , 



Figure lA.ll compares the electron capture rate on ^^Fe computed using the 
shell-model calculations of [22] with and without screening corrections. In ad- 
dition, the approximation (IA.13I) is also shown. The difference between the 
exact implementation of screening corrections, eq. flA.lOp . and the approxi- 
mation is never larger than 20% for the conditions shown in figure lA.ll 
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